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Electromagnetic wave scattering from a perfectly reflecting self-affme surface is considered. 
Within the framework of the Kirchhoff approximation, we show that the scattering cross section can 
be exactly written as a function of the scattering angle via a centered symmetric Levy distribution 
for general roughness amplitude, Hurst exponent and wavelength of the incident wave. The ampli- 
tude of the specular peak, its width and its position are discussed as well as the power law decrease 
(with scattering angle) of the scattering cross section. 
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Wave scattering from rough surfaces has been studied 
for long plpl (see also ||-[|] for recent reviews of the sub- 
ject) with potential applications in remote sensing, acous- 
tical/optical/radar detection or design of surfaces with 
specified scattering properties... The theoretical predic- 
tion of the angular distribution of the scattered intensity 
requires a proper statistical description of the surface 
roughness. At this stage it is classical to assume that 
the statistics of the height and its correlation function 
are Gaussian. This assumption allows to build analytical 
expressions and to compare them with experimental or 
numerical results. 

Since the pioneering work of Mandelbrot however, 
scale invariance has emerged as a relevant tool to de- 
scribe the geometry of real "disordered" objects. In the 
case of surfaces, the scale invariance takes the form of 
self- affinity. Such surfaces remain invariant under dila- 
tion of different ratios over the horizontal and the ver- 
tical directions. This long range correlated height dis- 
tribution is characterized by a roughness exponent, an 
amplitude parameter, and by the lower and upper limits 
of the scaling invariance region (instead of a height stan- 
dard deviation and a two point correlation function of 
finite width). It turned out that many real surfaces can 
be described through this formalism. Surfaces obtained 
from fracture growth or deposition processes |l(]] 

are classical examples. Similar results have been more re- 
cently found for surfaces obtained from cold metal rolling 
pd| . Note that such scaling invariance properties are ex- 
hibited from nanometric scales for native float glass |l^] 
up to kilometric scales for geological faults p3| . 

The aim of this letter is to exhibit an analytical expres- 
sion of the distribution of the scattering cross section for 
a perfectly reflecting self-affine surface in the framework 
of the Kirchhoff approximation. Since a first paper ex- 
amining the effect of scale invariance in the context of 
scattering from rough surfaces [ M , a large amount of 
studies have been published in various journals (see for 
example Ref . |l5|-|23| ) . Two main points motivated most 
of these studies i) the effect of scale invariance on the 



field scattered from the surface in comparison with the 
more classical case of a finite width correlation length ii ) 
the hope of directly measuring the "fractal" parameters 
of the surface via an acoustical or optical experiment. In 
spite of a large amount of works, it is fair to say that 
concerning far field scattering the only stable result that 
emerged is the existence of a power law tail in the in- 
tensity distribution with an exponent directly related to 
the roughness exponent. Very few analytical predictions 
can be found and oppositely a large number of numerical 
simulations still lack a clear physical interpretation. 

A reason for this could be an excessive focus on the 
"fractal dimension" (conversely the roughness exponent) 
with a blind point on the other parameters (the limits of 
the scale invariance regime and more often the amplitude 
parameter). A noteworthy exception is, however, due to 
Jakeman and his collaborators |24|,|2^] who predicted a 
Levy distribution for the intensity of a wave scattered 
from a self-affine random phase screen. However, these 
authors did not extend these analytical results to the case 
of scattering from surfaces. 

In the following we briefly recall the main properties 
of self-affine surfaces and we use the Kirchhoff approxi- 
mation to show that the scattering cross section can be 
written as a function of the scattering angle via a cen- 
tered symmetric Levy distribution. 

A surface is self-affine between the scales £_ and 
if it remains (cither exactly or statistically) invariant in 
this region under transformations of the form: 

(x,y,z) -> (Xx,Xy,X H z) . (1) 

The exponent, H, is usually called the roughness or 
Hurst exponent. Restricting ourselves to profiles, a sta- 
tistical translation of the previous statement is that the 
probability p(Az; Ax) of having a height difference Az 
over the distance Ax or its cumulative V(Az; Ax) = 
J- Z p($ z ', Ax)dSz is such that: 

V(Az;Ax)=V(X H Az-XAx) =${-£l\. (2) 
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Simple algebra based on this scaling gives 
a {Ax) =e 1 ~ H Ax H , 



(3) 



where the plane incident wave is given by: 
<&o{x, z\uj) = exp {ikx — iao(k, to)z} 



where a is the standard deviation of the height pro- 
file. Here I denotes a length scale, also known as the 
topothesy. This quantity is defined by o~{£) = I, which 
allows the geometrical interpretation of the topothesy as 
the length scale over which the profile has a mean slope 
of 45 degrees. The smaller £, the flatter the profile ap- 
pears on a macroscopic scale. In the case of a Gaussian 
height distribution, the probability p(Az; Ax) reads: 
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The self-affinc profile is thus fully characterized by its 
exponent H, its topothesy parameter I and the bounds 
of the self-affine regime £_ and £+. Let us mention in 
addition that the above geometric interpretation of the 
roughness amplitude made above does not imply that I 
has to lie between the lower (£_) and upper (£ + ) cut- 
off for the self-affine regime. For the surfaces usually 
considered in scattering problems, we rather expect that 
I <C When £_ < I < £ + , the topothesy makes the 
transition between the scales, below £, for which a frac- 
tal dimension D = 2 — H can be measured using the box 
counting method and the scales, above £, for which this 
dimension is just unity. 

In the following we consider the scattering of s- 
polarizcd electromagnetic waves from a one dimensional, 
random, self-affine surface. This surface, of Hurst ex- 
ponent H 7 we denote by z = C(x), and it is assumed 
to be Gaussian self-affine. We will further assume that 
the lower limit of the self-affine regime is smaller than 
the wavelength, A, of the incident wave. The scattering 
geometry considered is depicted in Fig. [I]. The incident 
plane is assumed to be the xz-plane, and the rough sur- 
face, which is perfectly conducting, is illuminated from 
the vacuum side by a plane wave of frequency u> = 2n/X. 
The incident and scattering angle respectively we denote 
by 6*o and 9, and they are defined positive according to 
the convention indicated in Fig. [|. In the above scatter- 
ing geometry there is no depolarization, and the electro- 
magnetic field is represented by the single non- vanishing 
component of the electric field $(x,z\u)) = E y (x,z\w), 
which should satisfy the (scalar) Helmholtz equation: 

'% + 8£ + ?z)*{x,z\w) = 0, (5) 



with vanishing boundary condition on z ~ and out- 
going wave condition at infinity. In the far field region, 
above the surface, the field can be represented as the sum 
of an incident wave and scattered waves: 



(7) 



and R(q\k) is the scattering amplitude. In the above 
expression, we have defined cto(q, u>) — i^/c) 2 — q 2 
(^tao(q,uj) > 0, ^sao(q,uj) > 0). Furthermore, the mo- 
mentum variables q and k are in the radiative region 
related to respectively the scattering and incident an- 
gle by q = (uj/c)sm9 and k = (u)/c) sin 9q, so that 
ao(q, u>) — {uj/c) cos9 and ao(fc, u>) — (co/c) cos6*o. 

Our aim is to obtain an expression for the scattering 
amplitude, since this quantity determines the scattered 
field. The mean differential reflection coefficient (or mean 
scattering cross section), which is an experimentally ac- 
cessible quantity, and defined as the fraction of the total, 
time-averaged, incident energy flux scattered into the an- 
gular interval (9,9 + d9), is related to this quantity by 
the following expression p6|: 
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Here L denotes the length of the self-affine profile func- 
tion as measured along the x-direction, and the other 
quantities have been defined earlier. The angle brack- 
ets denote an ensemble average over the surface profiles 
^(x), and the momentum variables are understood to be 
related to the angles 9q and 9 according to the expres- 
sions given above. 
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FIG. 1. Sketch of the scattering geometry. 

We now use the Kirchhoff-approximation, which con- 
sists of locally replacing the surface by its tangential 
plane, and thereafter using the (local) Fresnel reflection 
coefficient for the local incident angle to obtain the scat- 
tered field. Notice here that dealing with a surface whose 
scaling invariance range is bounded by a lower cut-off £_ 
does ensure that the tangential plane is well defined in 
every point. The scattering amplitude is then pqj: 



®{x,z\u>) = $ (x,z\u) 



—R(q\k) e^ x+iao ^ z , (6) R{q\k) = 
2ir 



2a (q,uj) J_ L/2 
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dx e'^-^^^Nix^), (9) 
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where j\T(x\lu) is a source function denned by J\f(x\u>) = 
2 d n $>o(x, ^\^>)\ z= ^f x y with the (unnormalized) normal 
derivative given by d n = —('(x)d x + d z . By substituting 
this expression for the source function into Eq. (||), one 
can show after some straightforward algebra where one 
takes advantage of the self-affine profile function having 
stationary increments, that the mean differential reflec- 
tion coefficient is given by 



exponent H. Figure g shows the mean differential reflec- 
tion coefficient (DRC), obtained for a incident angle of 
6 = 50° for different values of the ratio £/X. Note that 
the Levy distribution obtained above is likely not to be 
specific of the Kirchhoff approximation: the integrand in 
( |l0| ) is indeed nothing but the reflection of the statistical 
translational invariance of the problem. 
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with A£(v) = C(x) — C(x + v). Note that the statistical 
properties of the profile function, ((x), enters Eq. ( |To| ) 
only through fi(u). With Eq. (^), one may now calculate 
the ensemble average contained in Q(v). For a Gaussian 
self-affine surface one gets: 



/oo , 
dz exp < — i— ( 
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By performing the change of variable: 
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in Eq. (|l0|), and letting the length of the profile extent 
to infinity, L — ► oo, one finally obtains the following ex- 
pression for the mean differential reflection coefficient: 
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where a = v^fcos ^ cos , and C a {x) is the cen- 
tered symmetrical Levy stable distribution of exponent 
a defined as 
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Note that in the above expressions the wavelength, A, 
only comes into play through the combination (i/X) 1 ^^ 
which appears both in the prefactor and in the argument 
of the Levy distribution. This quantity can be geomet- 
rically regarded as s(A) 1/,ff where s(A) = {l/\) x ~ H is 
the typical slope of the surface over a wavelength. The 
behavior of the scattered intensity is thus entirely de- 
termined by this typical slope s(A) and the roughness 
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FIG. 2. The mean DRC as denned by Eq. (|Fj) for an inci- 
dence angle of 8q — 50° and a roughness exponent of H — 0.7. 
The surfaces had topothesy over wavelength of lj\ = 10 -3 , 
10" 4 , 10" 5 and 10" 6 which gives typical slopes at the wave- 
length scale of respectively s(A) = 0.126, 0.063, 0.032 and 
0.016. 



Taking advantage of the asymptotic expansion of the 
Levy distribution around zero [ p7| we find that the am- 
plitude of the specular peak scales as 
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and its half width (assuming — IcosOq <C 1) as 
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Notice that in the case of a non zero angle of incidence, 
9o, this (specular) peak is located at 8 = 8q + A8q, where 
A8q (A6>q ~ w 2 <c w) scales as 
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i.e. we have, due to the self-affinity of the surface, a 
shift, A8q, of the specular direction as compared to its 
"classical" position 8 = 8q. 

It is interesting that these non-trivial scaling results 
can all be retrieved via simple dimensional arguments. 
Let us examine the intensity scattered in direction 8; in 
a naive Huyghens framework two different effects will 
compete to destroy the coherence of two source points 
on the surface i) the angular distance separating 8 from 
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the specular direction ii) the roughness. The two source 
points will interfere coherently if the distance between 
them is lower than the two length scales d an g = A/p 
where p = tan[(0 — )/2] is the slope of the isophase 
plane and d roug h , the typical length scale over which the 
height difference is of order A. In case of a self-affine sur- 
face, d roug h = Xh^-h. Depending on the observation 
angle 9, the coherence will be controlled either by d roU gh 
or by dang, which will correspond respectively to the peak 
or the tail of the intensity distribution. The width w of 
the peak can then be defined as the transition between 
these two regions and will be such that d roug h — d ang . 
This leads directly to w oc (l/X)^^ 1 . While discussing 
the width of the peak present in the mean DRC, we also 
have to take into account the fact that real self-affine 
surfaces are not scale invariant over an infinite range of 
length scales. The upper cut-off £ + introduces a con- 
volution of the preceding results by a function of width 
A/£ + . This convolution does not affect the tails of the 
distribution but considering the peak we have now to dis- 
tinguish two regimes whose transition will be defined by 



(Contract No. 32690/213), Norsk Hydro ASA, and Total 
Norge ASA for financial support. 



X—HtH 



A where a denotes here the standard de- 



viation of the height profile measured over a macroscopic 
scale i.e. larger than £+. For a < A the width of the 
peak will be w ~ (•^/A)«~ 1 ~ s(\)~B and for <r > A we 
will get w ~ A/£ + . 

Using now the expansion of the Levy distribution at 
infinity |28| we get for the diffuse intensity (6 ^ do). 

^ T(l + 2H) sin(Trg) (M) 2 ~ 2g (cos °±^f- 2H 

47T - I • 



\d9 , 



cos o | sin <^| 1+2H 



The dependence of the prefactor on the surface parame- 
ters I and H is a new result of this letter. If the scaling 
can be easily retrieved via a first order expansion in plane 
waves, the above expression suggests that provided single 
scattering is considered, the power law behavior survives 
even in a non perturbative approach. 

We proposed a direct derivation of the intensity distri- 
bution of a wave scattered from a perfectly conducting 
self-affine rough surface in the framework of single scat- 
tering. Using a simple Kirchhoff approximation we ob- 
tain an exact solution consistent with all previously pub- 
lished perturbative or numerical results. Among other 
examples, our results could be directly tested in the case 
of an acoustic wave scattered from a fracture surface. 
This work will be extended to deal with more compli- 
cated problems directly related to scattering from rough 
surfaces (shadowing effects, multiple scattering, dielec- 
tric media...). Numerical investigations about the range 
of validity of the results obtained in this letter will follow. 

The authors would like to thank Claude Boccara, 
Jean-Jacques Greffet, Tamara Leskova and Alexei A. 
Maradudin for useful comments about this work. I.S. 
would like to thank the Research Council of Norway 



[1 
[2. 

[3 

[5; 

[6 

[-: 

[8 

[9 

[10 

[11 
[12 

[13 
[14 

[is; 

[16 

[17; 
[is; 

[19 

[20 
[21 

[22 

[23 

[24 

[25 

[26 

[27 

[28 



J.W.S. Rayleigh "The theory of sound", Dover, (New- 
York, 1945). 

P. Beckmann and A. Spizzichino, "The scattering from 
electromagnetic waves from rough surfaces", Artech 
House, (1963). 

J. A. Ogilvy, Theory of wave scattering from random 
rough surfaces", IOP Pub., (Bristol, GB, 1991). 
M. Nieto-Vesperinas, "Scattering and diffraction in phys- 
ical optics", Wiley, (1991). 

A. G. Voronovich, "Wave scattering from rough surfaces", 
Springer Verlag, (1994). 

B. B. Mandelbrot, "The fractal geometry of nature" , Free- 
man (1975). 

B. B. Mandelbrot, D.E. Passoja and A.J. Paullay, Nature 
308, 721, (1984). 

E. Bouchaud, E. Lapasset and J. Planes, Europhys. Lett. 
13, 73, (1990). 

K.J. Mal0y, A. Hansen, E.L. Hinrichsen and S. Roux, 
Phys. Rev. Lett. 68, 213, (1992). 

P. Meakin, "Fractals, scaling and growth far from equi- 
librium", Cambridge University Press, (1998). 

F. Plouraboue and M. Boehm, Trib. Int. 32, 45, (1999). 
D. Abriou, Y. Levillain, H. Arribart and F. Creuzet, 
Verre 2, 3, (1996). 

J. Schmittbuhl, Ph.D. thesis, Universite Paris 6, (1994). 
M. V. Berry, J. Phys. A 12, 781-797, (1979). 

D. L. Jaggard and X. Sun, J. Appl. Phys. 68, 11, 5456, 
(1990). 

M.K. Shepard, R.A. Brackett and R.E. Arvidson, J. Geo- 
phys. Res. 100, E6, 11709, (1995). 

P.E. McSharry, P.J. Cullen and D. Moroney, J. Appl. 
Phys. 78, 12, 6940, (1995). 

N. Lin, H.P. Lee, S.P. Lim and K.S. Lee, J. Mod. Optics 
42, 225, (1995). 

J. Chen, T.K.Y. Lo, H. Leung and J. Litva, IEEE Trans. 
Geosci. 34, 4, 966, (1996). 

C. J.R. Sheppard, Opt. Comm. 122, 178, (1996). 

J. A. Sanchez-Gil and J. V. Garcfa-Ramos, Waves Ran- 
dom Media 7, 285-293, (1997). 

J. A. Sanchez-Gil and J. V. Garcfa-Ramos, J. Chem. 
Phys. 108, 1,317, (1998). 

Y-P Zhao, C.F. Cheng, GC. Wang and T.M. Lu, Surface 
Science 409, L703, (1998). 

E. Jakeman, in "Fractals in physics", L. Pietronero and 
E. Tossati Ed., Elsevier, (1986). 

D. L. Jordan, R.C. Hollins, E. Jakeman and A. Prewett, 
Surface topography 1, (1988). 

A. A. Maradudin, T. Michel, AR. McGurn, and E.R. 
Mendez, Ann. Phys. 203, 255, (1990). 



At zero, £a(z) = ^r (±) 



At infinity C a (x) = 



r(i+a) 

tt\x\ 1 + " 



+ 0(x 4 



4 



